Universal renormalization-group dynamics at the onset of chaos in logistic maps and 

nonextensive statistical mechanics 
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We uncover the dynamics at the chaos threshold fioo of the logistic map and find it consists of 
trajectories made of intertwined power laws that reproduce the entire period-doubling cascade that 
occurs for \i < [1^. We corroborate this structure analytically via the Feigenbaum renormalization 
group (RG) transformation and find that the sensitivity to initial conditions has precisely the form 
of a g-exponential, of which we determine the q- index and the g-generalized Lyapunov coefficient X q . 
Our results are an unequivocal validation of the applicability of the non-extensive generalization of 
Boltzmann-Gibbs (BG) statistical mechanics to critical points of nonlinear maps. 
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Critical points of nonlinear maps offer a suitable play- 
ground for testing the validity of the non-extensive gen- 
eralization of the Botzmann-Gibbs (BG) statistical me- 
chanics proposed by Tsallis over a decade ago [1,2]. Here 
we describe universal properties related to the dynam- 
ics of iterates at the onset of chaos in unimodal maps [3] , 
that provide a literal confirmation of the generalized non- 
extensive theory. To this end we employ the celebrated 
one-dimensional logistic map, f^(x) = 1 — /i|s| 2 , —1 < 
x < 1, and the properties of its renormalization group 
(RG) fixed point, to present evidence of previously un- 
exposed scaling properties at the onset of chaos fi = fx^ . 
At this state, the most prominent of the map critical 
points, the trajectories of the iterates exhibit an intri- 
cate structure, that we describe and show is governed by 
the Feigenbaum's RG transformation [3]. 

The domain of validity of BG statistical mechanics has 
been implicitly challenged by the proposal of its non- 
extensive generalization. Subsequent studies have of- 
fered experimental and numerical evidence that point 
out both the inadequacy of the standard BG statistics 
and the plausible competence of the generalized theory 
in describing various types of phenomena and systems. 
This theoretical development represents an exceptional 
event in the long and trustworthy history of BG statisti- 
cal mechanics. However, it is still in the process of being 
converted into a rigorously corroborated and fully un- 
derstood fact. The suggested circumstances under which 
the generalized theory is believed to be applicable, at 
least with regards to non-linear dynamical systems, are 
those associated to a phase space with power-law sen- 
sitivity to initial conditions, to the consequent vanish- 



ing of the largest Lyapunov exponent, and to a fractal, 
or multifractal geometrical structure [2]. Here we show 
that our results for the dynamics at the onset of chaos 
in unimodal maps constitute an unequivocal proof of the 
universal validity of the non-extensive statistics at such 
critical points. 

In fact, at the chaos threshold (as well as at other 
critical points of the map) the Lyapunov exponent Ai 
vanishes, and the sensitivity to initial conditions £t, for 
large iteration time t, ceases to obey exponential behav- 
ior, exhibiting instead power-law behavior [4]. In order 
to describe the dynamics at such critical points, the q- 
exponential expression 



& = exp,(A,t) = [l-(q- l)A,t] 



-1/(9-1) 
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containing a ^-generalized Lyapunov coefficient X q , has 
been proposed [5]. This expression is based on the non- 
extensive entropy of Tsallis [2]. As a companion to this, 
generalizations for the Kolmogorov-Sinai (KS) entropy 
K q and for the Pesin identity A g = K g , X q > have 
also been introduced [5] (the standard expressions are 
recovered when q — > 1). Several recent studies [5-8], 
that probed numerically the onset of chaos of the lo- 
gistic map and its generalization to nonlinearity £ > 1 
(fn.({ x ) = 1 — A* \ x \ i — 1 < ^ < 1), have revealed a series 
of precise connections between the Tsallis entropic index 
q and the map basic parameters. Here we present RG 
analytical results that corroborate the previously known 
value of q at fioc for ( — 2 and also determine A g for the 
first time. 

To state our results more precisely, we recollect the fol- 
lowing properties. The logistic map exhibits several types 
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of infinite sets of critical points that appear as its control 
parameter /i varies; these correspond, amongst others, 
to period-doubling and chaotic-band-splitting transitions 
[3] . The accumulation point of the period doublings and 
also of the band splittings is the Feigenbaum attractor 
that marks the threshold between periodic and chaotic 
orbits, at fi^ — 1.40115... . The locations of period 
doublings (at /i = /i n < /.too) and band splittings (at 
/j = fi n > fioc) obey, for large n, power laws of the form 
Lin—ftcc ~ $~ n and ^00— M« ~ $~ n , where 5 — 4.6692... is 
one of the two Feigenbaum's universal constants. For our 
use below, we recall also the sequence of parameter values 
~p n employed to define the diameters of the bifurcation 
forks d n that form the period-doubling cascade sequence. 
At fi = ~p n the map displays a 'superstable' periodic orbit 
of length 2™ that contains the point x = 0. For large n, 
the distances to x = of the iterate positions in such 
2™-cycle that are closest to x = 0, d n = /A (0), have 
constant ratios: d n /d n +\ = —a, where a = 2.50290... is 
the second of the Feigenbaum's constants. A set of di- 
ameters with scaling properties similar to those of d n can 
also be defined for the band splitting sequence [3]. For 
clarity of presentation of our results we shall only use ab- 
solute values of positions, so that the dynamics of iterates 
do not carry information on the self-similar properties of 
"left" and "right" symbolic dynamic sequences [3]. This 
choice does not affect results on the sensitivity to initial 
conditions. Below, d n means \d n \. 



The main points in the following analysis are: 1) The 
iterates at fi^ follow trajectories that proceed in a con- 
certed manner according to the entire period-doubling 
cascade that takes place for fx < //qq. The positions of 
the trajectories are given in fact in terms of the diameters 
dn(~p n ) 01 the 2"-supercycles. 2) As a consequence, the 
sensitivity to initial conditions also evolves in agreement 
to the period-doubling cascade. 3) The bounds, or en- 
velops, as well as other monotonic subsequences, of both 
a single trajectory x t and of the sensitivity to initial con- 
ditions £t have precisely the form of a g-exponential. For 
£t we have q = 1— In 2/ In a and X q = In a/ In 2. 4) These 
results are obtainable via the fixed-point solution g(x) of 
the RG doubling transformation, consisting of functional 
composition and rescaling: TLf(x) = af(f(x/a). 
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FIG. 1. Absolute values of positions of the first 10 itera- 
tions r for two trajectories of the logistic map with initial 
conditions xo = (empty circles) and xo = S — 5 10~ 2 (full 
circles). Plotted quantities are dimensionless. 

To begin, we show in Fig. 1 the absolute values of the 
positions of two trajectories of the logistic map with ini- 
tial conditions xq = and xq = S ~ 5 10~ 2 , for the first 
10 iterations r. For r = 1 the positions are x\ = 1 and 
x\ ~ 1 — 3.5 10~ 3 , and it can be observed that the dif- 
ference between the two positions at times t = 2, 4, and 
8 grows progressively. In Fig. 2 we show the same first 
trajectory xq — and a second one xq = S ~ 10~ 4 , up 
to r = 1000. In the logarithmic scales it can be clearly 
appreciated that they consist of interwoven monotonic 
position subsequences with power-law decay. We are in- 
terested in the position subsequences that are generated 
by the time subsequences r = 2™ + 2™~ fe , k = 0, 1, ... . 
As we show below, the values for the trajectory subse- 
quence t = 2 n , with xq = 0, are asymptotically given 
by x T = d n , n > 0. More generally, each position sub- 
sequence t = 2" + 2"~ fc , obtained with xq = 0, is given 

by 5 (2k+1) (0) d k n , where <f£ = ^"^{Q) and n > k. 
The itinerary of the iterate starting at x$ = can be 
clearly observed in Fig. 2. The position approaches the 
origin xq — progressively as n increases every time that 
r = 2™, but in between the values 2™ and 2™ +1 it returns 
in an oscillatory manner towards x% = 1, repeating twice 
the positions visited in the previous cycle between 2" _1 
and 2™ and introducing a new position between these two 
sub-cycles. For small n the positions are approximately 
repeated, but they become accurately reproduced as n 
increases. The whole time series has the period-doubling 
structure. 
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FIG. 2. Absolute values of positions of the first 1000 iter- 
ations t for two trajectories of the logistic map with initial 
conditions xq = (empty circles) and xo — 5 — 10~ 4 (full 
circles) in logarithmic scales. The power-law decay of several 
time subsequences can be clearly appreciated. 

To interpret the dynamics in Figs. 1 and 2 in terms of 
the RG transformation, we consider R applied n times 
to the fixed-point map g(x), i.e.: 

g(x) =RWj(i) = a n gV n \x/a n ). (2) 

We determine the trajectory positions at times r = 2", 
with xo = 0. Since g(0) = 1, we have </ 2 ™)(0) = u~ n and 
because d n /d n +i — a with do — 1 implies d n — a~ n , we 
also have d n = g( 2 )(0). Thus, we obtain the diameters 
d n to be the positions X2™ • This result can be expressed 
as a g-exponential if we shift the time variable by one 
unit, t = 2™ - 1, and rearrange a~ n as (1 + t)- lna / ln2 . 
We obtain 

xt — expg(AQt), (3) 

with Q = l + ln2/lna and Aq = — In a/ In 2. Other po- 
sition subsequences t = 2™ + 2 n_fc can be put in the form 
of q-exponentials with the same values of Q and Aq. 

The expression for the sensitivity to initial conditions 
can be derived with the use of the following approximate 
property (that becomes asymptotically exact in the limit 
n — > oo) 

9 (2n) (dj) = — - =d n - ^d 2j - n , n < j. (4) 

To prove the above, note first that g^ 2 \dj) = 
g(2"+2 J )( ) = a~™ 5 ( 2 ^" +1 )(0), or, since g^" +1 \Q) = 
g(dj- n ), g( 2 >(dj) = a~ n g(dj- n ). The preceding equal- 
ity, together with g(dj- n ) = 1— fj,ood 2 _ n = 1 — 
Mooa 2 " -2 "', yields Eq. (4). Now, the distance between 



positions X2^(dj) — g( 2 '(dj) and x%n(di) = g^ n \di) at 
time r = 2™ can be written with the use of Eq. (4) as 
X2*i(dj) — X2"(di) — [x 2 a(dj) — X2o(di)]a n , or with use of 
the shifted time variable t — 2 n — 1 as 

x t {dj) - x t (di) = (x (dj) - x (di))a n . (5) 

The sensitivity to initial conditions £t is defined as 

£ t = lim \x t (dj) - x t (di)\ / \x (dj) - x (di)\ (6) 

|Aao|— >0 

(where limiAx i— >o is equivalent to limj ; j_ >00) j^y ) . £ t can 
then be written, considering that a n = (1 + t) lna ' ln2 , as 
the ^-exponential 

6 = exp g (A g i), (7) 

where g = 1 — In 2/ In a and X q = lncv/ln2. Notice that 
q = 2 — Q as exp g (y) = l/expg(— y). The previous con- 
struction applies strictly to initial positions that lie on 
the attractor; nevertheless, we remark that all other posi- 
tions tend to the attractor with a power-law behavior (see 
[10,11]), so that after an initial transient their dynamics 
becomes practically indistinguishable from the situation 
we describe here. In Fig. 3 we show the g-logarithm of £t 
vs t (with q = 1 — In 2/ In a = 0.2445...), from a numer- 
ical simulation of two trajectories with initial conditions 
xq = and xo — S ~ 10~ 8 . The result is a straight 
line with slope very close to X q = In a/ In 2 = 1.3236... . 
This corroborates the RG prediction (the g-logarithm, 
ln 9 y = (y 1 ^' 1 — !)/(! — q), is the inverse of exp (y)). 
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FIG. 3. The g-logarithm of sensitivity to initial conditions 
£t vs t, with q = 1 — ln2/lna = 0.2445... and initial con- 
ditions xq = and xo = S ~ 10~ 8 (circles). The full 
line is the linear regression y(t). As required, the numeri- 
cal results reproduce a straight line with a slope very close to 
X q = In a/ In 2 = 1.3236... . 



3 



Interestingly, both x t /xo and £ t can be seen to satisfy 
the dynamical fixed-point relations h(t) — ah{h{t/a)) 
with a = 2 1 /W~ 1 ) and a = 2 l ^ c '- 1 \ respectively. In 
relation to this, we note that the static fixed-point so- 
lution f*(x)/x to the Feigenbaum RG recursion relation 
for the case of the tangent bifurcation, obtained by Hu 
and Rudnick [3,9] to study the intermittency transition in 
the ^-logistic map, has the form of a g-exponential with 
q = 2. This has been pointed out recently [10,11], where 
in addition it has been shown that this solution applies 
too, but now with q — 3, to the period-doubling transi- 
tions that take place at fi n < Moo- We recall that, for the 
transition to periodicity of order n, the RG transforma- 
tion is applied to the n-th composition /(") of the original 
map in the neighborhood of one of the n points tangent 
to the line with unit slope, and a shift is made of the 
origin of coordinates to that point. The RG fixed-point 
map is f*(x) = xexp^ua; 9-1 ), where u is the leading 
expansion coefficient of and where the recursion re- 
lation is satisfied with a = 2 1 ^ q ^ 1 \ In the neighborhood 
of the intermittency and period-doubling transitions the 
time evolution of iterates follow monotonic paths set by 
the form of the map itself, and are also g-exponentials. 
For these types of critical points the static and dynamic 
properties are simply related and obey expressions of the 
same form [10,11]. Not only the entropic index q can be 
plainly identified, but the ^-generalized Lyapunov coeffi- 
cient X q turns out to be given by the expansion coefficient 
u. Unequivocal corroboration of these results has been 
obtained recently [11]. 

At the onset of chaos both static and dynamical prop- 
erties are more complex. The celebrated RG fixed-point 
map static solution is obtained as a power series of a 
smooth unimodal trascendental function, the universal 
n — > oo limit of (— a) n /l 2 ^ (x/(— a) n ); likewise, the frac- 
tal dimension of the attractor df — 0.538045143... is ob- 
tained considering also the same n — + oo limit on posi- 
tions of the 2™-cycles [3,12]. As we have seen here, the 
multifractal attractor at /Hoc imprints an involved struc- 
ture into the time evolution of the iterates, that can be re- 
solved in terms of simpler monotonic time subsequences. 
Remarkably, these subsequences and the sensitivity to 
initial conditions £ t are analytically reproduced by the 
same RG transformation originally applied to describe 
static properties. These quantities evolve as universal q- 
exponentials, with q and X q simply expressed in terms of 
a. We observe then a connection between dynamic prop- 
erties of a strange attractor at the edge of chaos, such 
as the ^-generalized Lyapunov coefficient, and the static 
properties, described by the set of distances d\ that make 
up this multifractal; this is in the spririt of the Kaplan- 
Yorke conjecture [3]. 

Under some conditions, exemplified here by critical 
points in nonlinear maps, the Lyapunov exponents of a 
system that measure the strength of phase space mixing 



vanish. When this happens dynamic processes become 
sluggish in exploring their permissible configurations and 
may be capable of covering only a small fraction of the 
available phase space, even in the limit t — > oo. This frac- 
tion may have a fractal dimension smaller than the total 
dimension of phase space. These are thought to be the 
conditions for failure of BG statistics and applicability of 
its nonextensive generalization [2] , and the significance of 
our analytically-backed results with no approximations is 
a contribution towards the clarification of this issue. At 
the onset of chaos of unimodal maps the reduced phase 
subspace is represented by the strange attractor, a Can- 
tor subset of the interval — 1 < x < 1. It is important 
to point out that in this case the permissible positions 
(configurations) are asymptotically confined by the at- 
tractor and this acts as an inescapable barrier to move- 
ment to other locations. By construction, the dynamics 
at the onset of chaos, as well as those restricted to the 
neighborhood of the intermittency transitions, describe a 
purely nonextensive regime. It should be mentioned that 
our analysis does not consider the access of trajectories 
to an adjacent or neighboring chaotic region, as in the 
setting of Refs. [13,14] or that in conservative maps [15]. 
Hence there is no feedback vehicle for a crossover from a 
nonextensive regime with vanishing ordinary Lyapunov 
exponent to an extensive regime with a positive one at 
some t cross as t — > oo. 

Clearly, our findings have wide-ranging validity, as 
they apply to all dissipative systems in the Feigenbaum 
universality class. It might be possible to observe anal- 
ogous behavior at chaos boundary criticalities in other 
classes of dynamical systems, such as period doubling in 
bimodal [16] and other multiparameter maps [17], and 
quasiperiodicity in dissipative systems [18]. More gen- 
erally, the properties discussed here are likely to hold 
strictly for other types of systems or situations that pos- 
sess equivalent phase space limitations. Those systems 
for which experimental and numerical evidence has accu- 
mulated on BG statistics inadequacy and nonextensive 
statistics competency [2] warrant examination. The more 
that is learned on mechanisms and circumstances leading 
to a hindered phase space the clearer the physical under- 
standing of the applicability of the nonextensive theory 
will become. 
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